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Abstract. We provide a probabilistic analysis of the upwind scheme for d-dimensional 
transport equations. We associate a Markov chain with the numerical scheme and then 
obtain a backward representation formula of Kolmogorov type for the numerical solution. 
We then understand that the error induced by the scheme is governed by the fluctuations 
of the Markov chain around the characteristics of the flow. We show, in various situations, 
that the fluctuations are of diffusive type. As a by-product, we recover recent results due 
to Merlet and Vovelle [13] and Merlet [T2]: we prove that the scheme is of order 1/2 in 
L°°([0,T],L 1 (M. d )) for an initial datum in BV{R d ) and of order 1/2 - e, for all e > 0, in 
L°°([0,T] x R d ) for an initial datum in W 1 '°°(M. d ). Our analysis provides a new interpreta- 
tion of the numerical diffusion phenomenon. 

Resume. Nous proposons une analyse probabiliste du schema upwind pour les equations de 
transport en dimension d quelconque. Pour cela, nous associons au schema une chaine de 
Markov qui nous permet d'obtenir une formule de representation de type Kolmogorov pour 
la solution numerique. Nous comprenons alors que l'erreur due au schema est gouvernee 
par les fluctuations de la chaine de Markov autour des caracteristiques du transport. Nous 
montrons, dans des situations diverses, que ces fluctuations sont de type diffusif. Commc 
consequence, nous retrouvons des resultats recents de Merlet et Vovelle [13] et Merlet [T2] : 
nous montrons que le schema upwind est d'ordre 1/2 dans L°°([0, T], L 1 (R d )) pour une 
donnee initiale dans BV(R d ), et d'ordre 1/2 — e pour tout e > pour une donnee initiale 
dans W 1 ' ca (R d ). Cette analyse donne une interpretation nouvelle du phenomene de diffu- 
sion numerique. 
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1. Introduction 

This paper provides a new analysis of the upwind scheme for the transport problem in 
dimension d EN \ {0} 

, . f d t u(t,x) + (a(x),Vu(t,x)) = 0, (t,x) e [0,T] x R d , 

^ } 1 u(0,x) = u°(x), x G R d . 
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We assume a to be Lipschitz continuous, so that (11. ip is well-posed. Several different reg- 
ularity assumptions are made on u° in the following, among which u° G W 1 '°°(M d ) and 
u° G BV(M d ). In any case, the unique solution to (II. ip is u(t,x) = u°(Z(t,x)) where Z is 
the backward characteristic, i.e. the solution of 

f t Z(t,x) = -o(x), (t,z) G [0,T] x R d , 
1 J \ Z(0,x) = x, x G M d . 

The upwind scheme is a standard method to solve this problem in an approximate way 
(see for instance [6j). It is derived and described in Section [2] below. 

In dimension 1, the scheme is known to be first order consistent (with respect to the 
maximal cell diameter h) with the transport equation. Thus, it is first order convergent, for 
any u° G C 2 (lR d ), provided that a Courant-Friedrichs-Lewy (CFL) stability condition holds. 
For non-smooth initial data, the upwind scheme is just of order 1/2. This loss of convergence 
order is traditionally attributed to the dissipative character of the scheme. Up to now, the 
use of the word "dissipative" has been justified by the following fact: on a uniform mesh, 
the scheme is second order consistent with an advection-diffusion equation, the diffusion 
coefficient being first order with respect to h; as a consequence, the numerical error at time 
t is proportional to Vth for a non-smooth initial datum. 

In this paper, we provide another explanation of the diffusive behavior, which is valid on 
any general mesh in dimension d. We here interpret the numerical diffusion by means of a 
stochastic process. Let us briefly describe the basic idea. The approximate value given by 
the upwind scheme in a cell K at time step n + 1 is a convex combination of the approximate 
values at time step n in the cells neighboring K: 

LeT 

where T is the set of cells, and Pk,l G [0,1] with J2lgtPk,l = 1 (see Section [2] for the 
complete definition of the scheme). This convex combination allows a probabilistic interpre- 
tation: we can define a random sequence of cells {K n ) neN as a Markov chain with probability 
transition, from K to L, Pk,l- In this framework, the upwind scheme appears as the expec- 
tation of a random scheme associated with the chain (K n ) n >o. Precisely, the value u\ is the 
expectation of the value of u° in the cell K n when K is chosen as the starting cell of the 
chain. In a probabilistic way, we write: 

u\ = E K «J , 

the symbol K in the notation E# meaning that K = K. (See Theorems 13.11 and 14.11 ) In 
the theory of stochastic processes, the above identity is a backward Kolmogorov formula: it 
is the analogue of the representation formula of the heat equation by the Brownian motion. 
We then understand the chain (i^ n ) n >o as a random backward characteristic. 

Our main idea consists in analyzing the behavior of the random characteristic according 
to the following program. The first point is to show that the mean trend of the random 
characteristic coincides with the exact characteristic Z, solution to (II. 2p . The next step is 
to understand that the error of the numerical scheme is governed by the fluctuations of the 
random characteristic around the exact one. Heuristically, the order of the fluctuations is 
given by the central limit theorem: therefore, we expect them to be controlled, in a suitable 
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sense, by Ch 1 / 2 where C only depends on the datum a and the time t. The final step is to 
derive the 1/2 order of the scheme. 

Applying this program, we establish the 1/2 order in L°°([0,T], L 1 (R d )) for u° G BV(R d ). 
(See Theorem 15.101 ) For u° G W 1 '°°(M d ), we also prove that the scheme is of order 1/2 — e 
in L°°([0, T], L°°(M d )) for all e > 0. (See Theorem EE) In this last case, it is clear that 1/2 
is an upper bound (see [9], and [19] for the non-linear case), but the exact convergence order 
remains unknown. The reason why our estimate is better in the L l -'m. space norm may be 
explained as follows. Estimating the error in L l amounts to average the initial cell of the 
random characteristic. This additional averaging reduces the weight of the trajectories of 
the chain that are away from Z. 

Since the pioneering article of Kuznetsov [TT], in which the 1/2 order is established in 
the Cartesian framework (for linear an non-linear scalar equations), many papers have dealt 
with the rate of convergence of the upwind scheme. Let us briefly review them. 

For general scalar equations with a datum u° G BV(M. d ), Cockburn, Coquel and Le 
Floch [I], Vila [21], [2] and Chainais-Hillairet [3] prove the (non-optimal) 1/4 order in 
L°°([0, T], L 1 (IR d )) under slightly different hypotheses (and for several schemes, including 
the upwind one). 

For hyperbolic Friedrichs systems, Vila and Villedieu [22] derive a 1/2 order estimate in 
the L 2 ([0,T] x R d ) norm for if 1 (M d ) initial data. 

In the frame of the linear transport equation, which we are involved in, Despres [3] proves 
a 1/2 order estimate in the L°°([0, T], L 2 (R d )) norm in the case of H 2 (R d ) data. His proof 
relies on a precise study of the consistency of the scheme after several time steps. (It is 
indeed known that the scheme is not consistent at each time step on a general mesh.) For 
C 2 initial data, Bouche, Ghidaglia and Pascal [lj show the order 1 in the L°° norm, under a 
condition on the mesh that is related to consistency. At last, in recent works, 

• for an initial datum in BV(M. d ), Merlet and Vovelle [13] show the optimal estimate 
of order 1/2 in the L°°([0, T), L 1 ^)) norm, 

• for a Lipschitz continuous initial datum, Merlet [T2] shows the order 1/2 — e, for any 
e > 0, in the L°° norm. 

It is thus understood that our paper provides a new proof of the results obtained in [13] 
and [12J. Actually, our framework is slightly different since we do not assume the velocity 
a to be divergence- free, as done therein, but we assume it to be independent of time. We 
think that this does not make fundamental differences. Despite the similarity of our results, 
we insist on the fact that the arguments here are completely different. As said above, our 
proofs rely on the analysis of the stochastic characteristic (K n ) n > (that shall mimic the exact 
characteristic Z). In particular, we do not use energy estimates. (Except those of Despres 
and Bouche et al. based on the consistency of the scheme, all the papers mentionned above 
are built on energy or entropy estimates). 

Our paper is organized as follows. In Section [21 we state the framework of our analysis. 
In Section [31 we focus on the one-dimensional case to introduce, with great care, the notion 
of stochastic characteristic. By the way, we establish a refined estimate of the order of the 
scheme in the specific case where the velocity is constant and the mesh is regular. (See 
Proposition 13.51) . In Section @J we extend the probabilistic interpretation of the upwind 
scheme to the higher dimensional setting. We then provide a direct proof of the 1/2 order 
in L°°([0, T], L 1 (IR d )) in the following simple case: u° is assumed to be periodic, as well as 
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the mesh, and Lipschitz continuous. This section is the heart of the paper. Refining the 
strategy, we finally obtain in Section [5] the announced results. This last part is a bit more 
technical and relies on a concentration inequality for martingales, which is given in Annex, 
see Sectional 

2. Framework and useful notations 

Let {K} K( - T , the mesh, be a set of closed polygonal subsets of M. d with non-empty disjoint 
interiors such that M. d = [j Ke7 - K. The volume (<i-Lebesgue measure) of a given cell K G T 
is denoted by \K\. The supremum of the diameters of all the cells is denoted by h, i.e. 
h = sup^ gr diam(if). Two cells K and L are said adjacent if they aren't disjoint but have 
disjoint interiors. In this case, we write K ~ L. We assume that, for all pairs (K, L) 
of adjacent cells, the intersection K fl L is included in a hyperplane of M. d . The surface 
((d — 1)-Lebesgue measure) of the face K fl L is then denoted by \K H L\. 

Let At > be the time step of the method. The value u\ intends to approximate the 
mean value of u(nAt, •) in the cell K. The upwind scheme provides a way to compute 
such u K . It is easily obtained by integrating the divergence form of the PDE in (11. ip . 
d t u + div(cra) — wdiv(a) = 0, over [nAt, (n + 1) At] x K. We get 

(2.1) / u((n + l)At,x)dx - I u(nAt,x)dx 

J K J K 

p r(n+l)At p p(n+l)At 

+ / / (a(x),nK,L,) u(t,x)dtdx — / / u(t, x)div (a) (x)dtdx = 0, 

L~K J KnL J nAt J K J nAt 

where tlk,l is the unit normal vector on K fl L outward from K. From a numerical point 
of view, it then seems natural to compute both an approximate value u\ of the mean of 
uinAt, ■) in cell the K, i.e. 



u K w 



\K\ 1 / u{nAt,x)dx, 

J K 



and an approximate value L of the mean of u on the edge K fl L between the time steps 
n and n + 1, i.e. 

•(n+l)At 



Ar^KDLl- 1 / / u(t,x)dtdx. 

JKnL J nAt 



The quantity u\ L is called the numerical flux. Defining clk,l as the mean value of a on the 
edge K fl L, i.e. 

a,K,L = \K R / a{x)dx, 
Jkhl 

we get the following approximate version of (12. ip . 
v n+1 — v n ^ 

1^1 K K + Yl ( a ^, n K , L ) \K n L\u n KtL -u K J2 < a ^' \ K ^L\= 0. 

The upwind scheme considers the numerical fluxes u\ L as upwinded: u\ L = u\ for L e K + 
and L = u\ for L e K _ with 

K+ = {L~ K, (a KtL , n K>L ) > 0} , 
K- = {L~ K, {a K)L , n K)L ) < 0} . 
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This finally gives 

_ 71+1 ,, 

(2.2) \K\ K K + ^ (a Xii ,nx,L)|^nL|K-<) =0, (^)GNxT. 



At 

LeK- 



The numerical initial condition is usually taken as u° K = \K\ 1 f K u°(x) dx. It is straightfor- 
ward that the scheme satisfies the maximum principle under the condition 

_ ^ {a KtL) n KtL )\KC\L\ <x KeT 



L&K- 



This condition is called the Courant-Friedrichs-Lewy (CFL for short) condition and is as- 
sumed to be satisfied in all the paper. 

3. Analysis in Dimension 1 

For pedagogical reasons, we first investigate the one-dimensional framework. As an- 
nounced in Introduction, the velocity field a : R — > R is assumed to be bounded and to 
be /t-Lipschitz continuous. In particular, for any starting point x G R, the characteristic 
equation starting from x 

(3.1) d t Z(t, x) = -a(Z(t, x)), t> 0, Z(0, x) = x, 

admits a unique solution. Denoting by u° the initial condition of the transport equation, 
which is assumed to be /t-Lipschitz continuous in the whole section, the solution of the 
transport equation rewrites 

(3.2) u(t,x) = u°(Z(t,x)), (t,x) 6 1 + x K. 

In this section devoted to dimension 1, for every cell K e T, the volume (length) of K 
is denoted Axk- The edge value of a is defined as o,k,l = a{K H L). The constants "C" 
and "c" used below only depend on Halloo and k. They are always independent of At, of 
h = sup KeT Ax K , of the time index n and of the random outcome uj. In particular, the 
notation 0(x), for a given variable x, denotes a quantity bounded by Cx for some constant 
C only depending on Halloo and k. 

3.1. Probabilistic Interpretation. In the one-dimensional framework, the scheme has the 
form 

u° K = -J— [ u°(x)dx, K eT, 
f3 V K Jk 

1 j n+i V- a K rAt n , , , ^ a K rAt, 



LeK- K LeK- K 

and the following CFL condition is assumed to be in force 

Ik,lAi 
Ax K 



LeK- LeK- 
sumed tc 

(3.4) Ker. 

LeK~ 

The geometry of the mesh is simple: each cell K has two neighbors. When the velocity field 
a is non-zero in the cell K, there is one and only one cell L in K~ . If a is positive in K, it 
is the left one; of course, if a is negative, it is the right one. 
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We focus for a while on a given cell K. By the CFL condition (13. 4p . all the coefficients 

Pk,l = 7 for L E K , 

Ax K 

a K)L At 



Axk 



Pk.k = 1 + ^2 

L£K~ 

p K)L = for L E T \ (K- U K) , 

are non-negative and may be seen as probability weights. Henceforth, for a given time step 
n > 0, the right-hand side in (13.31) may be interpreted as an expectation with respect to 
these weights: 

Intuitively, this means that we are choosing one cell among K U {L ~ K}, K being fixed, 
with the probability weights pk,k and Pk,l for L ~ K. To make this idea rigorous, we 
introduce a probability space (Q, A, P) as well as a random variable £ : Q — > K U {L ~ X} 
such that P{£ = L} = p^.L for any L ~ K and P{£ = K} = Pk,k- Then, the (n + step 
of the numerical scheme on the cell K can be written in the following way: 

(3.5) u n K +1 = E>^< = E K1" 

This relationship provides a probabilistic interpretation for the one step dynamics of the 
numerical scheme. We are to iterate this procedure. 

The probabilistic dynamics between times n and n + 1 just depend on the starting cell 
K. In the theory of stochastic processes, this property is typical of Markovian dynamics. 
Indeed, the family of probability weights (pk,l)k,lgt defines a stochastic matrix of infinite 
dimension (all the entries of the matrix are non-negative and the sums of the entries of a 
same line are all equal to 1). This stochastic matrix corresponds to the transition matrix 
of a Markov chain. Up to a modification of the underlying probability space, there exists a 
sequence (i^ n ) n >o of random variables taking values into the set of cells as well as a collection 
of probability measures (f^KeT, indexed by the cells, such that, under each P#, (K n ) n > 
is a Markov chain with rates {pk,l)k,lgT starting from Kq = K. In other words, 

Vn > 0, F K {K n+1 = L\K n = K} = p KiL , F K {K = K} = 1. 

The behavior of the chain (i^ n ) n >o is as follows: if the velocity is positive in the cell K n , 
then the probability Pk„,l vanishes if L is the right neighbor of K n , so that the chain can 
either stay in K n or jump to the left. 

Now, we can interpret (13.51) in a different way: 



u 



n+l 
K 



where E# denotes the expectation associated with P#. This means that is the expec- 
tation of u n in the random cell K\ occupied by the Markov chain, which started one time 
step before in K. We can also write for any integer % > 



u 



n+l 



Kk [ux i+1 \Ko, . . . , Ki\ Px— almost surely. 



(When conditioning with respect to Kq, . . . ,Ki, the past before i — 1 doesn't play any role, 
and the chain restarts, afresh, at time i from Ki.) In what follows, we denote the conditional 
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expectation [-\Ko, ■ ■ ■ , K n ] by We also omit to specify that such a conditional 

expectation is computed under P#. Now, we are able to iterate the procedure in (13.51) : 

u n K +1 = E K [u n Ki ] =E K [E 1 K [u n K2 1 ]} 

We have proved the following representation for the numerical solution u n : 

Theorem 3.1. Under the above notations, the numerical solution u\ at time n and in the 
cell K has the form: 

u n K = E K [u° Kn ]. 

The representation given by Theorem 13.11 is a backward Kolmogorov formula for the nu- 
merical scheme. Generally speaking, the backward Kolmogorov formula provides a represen- 
tation for the solution of the heat equation in terms of the mean value of the initial condition 
computed with respect to the paths of the Brownian motion (or of a diffusion process). In 
this framework, the paths of the Brownian motion appear as random characteristics. In our 
own setting, the Markov chain (K n ) n > Q almost plays the same role. 

We say "almost plays" because the sequence (K n ) n > is n °t a sequence of points as the 
Brownian motion is. Actually, we have to associate with each random cell K n a random 
point X n (X n being ideally in K n ) to obtain a random characteristic (X n ) n > . 

The choice of these points is crucial. In what follows, we choose X n as the entering point 
in the cell K n . This means that 

X n = X n _i if K n = K n _i, X n = K n f] K n _i if K n ^ K n _ x . 

The above definition holds for n > 1. The position of the initial point X inside K has to be 
specified. If a(x) > for all x G K , we choose X as the right boundary of K . (Indeed, the 
right boundary plays in this case the role of the entering point since the velocity is positive.) 
If a(x) < for all x G K , we choose X as the left boundary. If 3x G K such that a{x) = 0, 
we choose X as the middle of K . 

What is important is that the sequence (X n ) n > is adapted to the filtration generated by 
(K n ) n >o, i.e. (a(Ko, . . . , K n )) n > : knowing the paths (K , . . . , K n ), one knows the positions 
of the points (X , . . . , X n ). 

The reader may wonder about this specific choice for the sequence (X n ) n > . Assume that 
the velocity a is non-zero, say for example positive, in the cell K n . By continuity, it is 
positive in the neighborhood of K n : the chain goes from the right to the left in this area of 
the space. As a by-product, the entering point in the cell K n is the right boundary of K n . 
In this case, A n+1 is either X n itself or the left boundary of K n , which is the right boundary 
of K~ (this set of cells is in the present case a singleton and we identify it with its element, 
as well as for in the following when the velocity is away from 0), so that 

E n K [X n+1 -X n ] = -Ax Kn p KnK - = -a KnK -kt. 

(Indeed, the probability that X n+ i is the right boundary of K~ is given by the probability of 
jumping from K n to K~.) In other words, the mean displacement from X n to X n+ i, knowing 
the past, is exactly driven by the velocity field —a. This is very important: loosely speaking, 
the speed of the random characteristic is given by the velocity field of the characteristic 
equation itself! Here is a more precise statement. 
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Proposition 3.2. For every n > 0, 

E n K [X n+l - X n ] = -a(X n )At + O(hAt). 

Proof. If a(x) > for all x G K n , then X n has to be the right boundary of K n . (By 
definition of X Q , this is true until the first jump of the chain. After the first jump, this is 
still true since the chain cannot come from the left by positivity of a.) Moreover, starting 
from K n , the chain cannot move to the right since a Kn K + = a(X n ) > 0. Hence, X n+1 has 
to be either X n or the left boundary of K n . As done above, 

E n K [X n+1 - X n ] = -Ax Kn p Kn>K - = -a KnjK -At = -a(X n )At + O(hAt) 

by the Lipschitz property of a. The same argument holds when a(x) < for all x G K n . 

If a(x) = for some x G K n , then a(X n ) = O(hAt) by the Lipschitz property of a. More- 
over, the probability of moving to the right is equal to max(a^ n K +, 0) x At/AxK n = O(At). 
The same holds for the probability of moving to the left. When moving, the displacement is 
bounded by h so that the result is still true. □ 

Remark. The necessity of choosing the point X n as the entering point in the cell K„ is 
related to the well-known fact that the upwind scheme is consistent (in the finite difference 
sense) with the transport equation provided that the control points for every cell are chosen 
on the right if the velocity is positive and, conversely, on the left if the velocity is negative: 
see [6]. 

3.2. A First Example: a and Ax constant. To explain our strategy, we first focus on the 
very simple case where both a and Ax are constant: a(x) = a for all x G K and Axk = h for 
all K G T '. Without loss of generality, we can assume that a is positive so that the random 
characteristic goes from the right to the left. In this setting, the transition probabilities are 
of the form 

VK G T, Vk,k- = —r-, Pk,k = 1 7— ■ 

h h 

The probability of jumping from one cell to another doesn't depend on the current state of 
the random walk. From a probabilistic point of view, this amounts to say that the sequence 
(X n+ i — X„) n > (with X equal to the right boundary of the initial cell) is a sequence of 
Independent and Identically Distributed (IID in short) random variables under Fx, whatever 
K is. The common law of these variables is given by 

W K {X n+1 -X n = -h} = l- P{X n+1 - X n = 0} = 

n 

In particular, we recover a stronger version of Proposition 13.21 ( "stronger" means that there 
is no 0(hAt)): 

E K [X n+1 - X n ] = -a At. 

In particular, the mean trend of the random characteristic is exactly driven by the velocity 
—a, that is by the mapping t \— > Xq — at, which corresponds to the characteristic of the 
transport equation with X as initial condition, i.e. Z(t,X ) (see (13. lft ). At this stage, we 
understand that the order of the numerical scheme is deeply related to the fluctuations of the 
random characteristic around its mean trend, that is around the deterministic characteristic. 
Indeed, for any starting cell K, we have 



u 



K 



E K [u° Kn ] = E K [u (X n )]+O(h), 
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where 0(h) only depends on the Lipschitz constant of the initial condition u° and is inde- 
pendent of the initial cell K. Thus 

n-1 

(3.6) u n K = E K [u°(X - anAt + ^(X fc+1 - X k + aAt))] + 0(h). 

k=0 

By (13T2|1 . for all x G K, 

n-1 

- u(nAt,x) = E K [u°(X - anAt + ^(X k+1 -X k + aAt))] - u(nAt,x) + 0(h) 

k=0 
n-1 

= E K [u°(X - anAt + ^(X k+1 - X k + aAt)) - u(nAt, X )] + 0(h) 

k=0 
n-1 

= E K [u°(X - anAt + ^(X k+1 - X k + aAt)) - u°(X - anAt)] +0(h). 

k=0 

Using again the Lipschitz continuity of u°, we deduce by the Cauchy-Schwarz inequality 
that, 

n-1 

\u n K -u(nAt,x)\ < nE K [\^2(X k+1 -X k + aAt)\] + 0(h) 

k=0 

n-1 

< KE K [\J2(X k+1 -X k + aAt) | 2 ] 1/2 + 0(h) 

k=0 

for all x £ K. The last expectation is nothing but the variance of the sum of the random 
variables (X k+X - X fe ) <fc<„.-i under F K , i.e. 

n— 1 n—1 

E K [\J2( X k+i ~ X k + «At)| 2 ] = ¥ K [J2( X k+i ~ X k )] ■ 

k=0 k=0 

It is well-known that the variance of the sum of independent random variables is equal to 
the sum of the variances of the variables. We deduce that, for all x e K, 

\u n K -u(nAt,x)\ < K[nV K (X 1 - X )] 1/2 + 0(h). 

The common variance is equal to 

Vjf(Xi - X Q ) = h 2 p K , K - - a 2 At 2 = aAt (a - aAt). 

Note that the CFL condition guarantees that the right-hand side above is non-negative, 
which ensures that the equality is meaningful. We thus recover a well-known estimate for 
the L°°-error induced by the upwind scheme: 

Proposition 3.3. Assume that a(x) = a and Axk are constant and that u° is bounded and 
K-Lipschitz continuous. Then, at any time n > 0, 

1 /2 

sup \\ux — u(nAt, ^Hioofjn — n(naAt(h — aAt)) + 0(h). 
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3.3. General Case: a(x) and Axk non constant. Our strategy is sharp enough to obtain 
the analogue of Proposition ^. 31 when a does depend on x and Axk on K. The main difference 
here is that Xq can be either the right boundary of Kq or the left boundary or the barycenter 
of K , according to the sign of the velocity in K . Following the previous subsection, the 
difference between the numerical and the true solutions at time n > on a cell K is given 
by 

u\ - u(nAt, x) = E K [u° Kn ] - u°(Z(nAt, x)) = E K [u° Kn - u°(Z(nAt, x))] 
for all x G K. As above, the Lipschitz property yields 

u° Kn =u°(X n ) + 0(h). 

Again, the term 0(h) is uniform with respect to the starting cell K, to the time index n, 
to the parameters At and h and to the underlying draw uj G Q. By Gronwall's lemma, we 
control the distance between Z(nAt,x) and Z(nAt, Xq), so that 

\u n K - u(nAt, x) | < kE k [\X n - Z(nAt, x) \] + 0(h) 

= nE K [\X n - Z(nAt,X )\] +0(h)exp(nnAt). 

We have 

(3.8) X n - Z(nAt, X )=X n -X + J2 <Z(s, X ))ds. 

k=0 JkAt 



(3.7) 



n— 1 n— 1 



k=0 k=0 



n—1 



By Proposition 13. 2\ 

(3.9) X n - X = J2( X ^ ~ X *) = ~&tJ2a(X k ) + M n + 0(nhAt), 
with 

(3.10) M n = ^(A fe+1 -X k - E k K (X k+1 - X k )) (M = 0). 

k=0 

By the boundedness and the Lipschitz continuity of a, 

/•(fc+l)At n-l 

(3.11) a(Z(s,X ))ds = AtJ2 a ( z ( kAt i X o)) +0(nAt 2 ). 
k=o J kAt k=Q 

Plugging <^M) and fl3~TT|) into (ET5I) . we obtain 

n-l 

|A„ - Z(nAt, A )| < «AiV |A fc - Z(fcAt,X )| + |M n | + 0(nhAt + nAt 2 ). 

k=0 

Taking the expectation of each term and applying Gronwall's lemma, 

(3.12) E K [\X n - Z(nAt,X )\] < [E K [\M n \] + 0(nhAt + nAt 2 )] exp(nnAt). 

As in the case where both the velocity and the spatial step are constant, the process 
(M n ) n > represents the fluctuations of the random characteristic around a discretized version 
of the deterministic characteristic. (See f)3.9p .) In the probabilistic theory, it is a martingale 
on (Q, A, ^k), i-e., at any time n > 0, M n is a(K , . . . , i^ n )-measurable and E^[M n+1 ] = M n . 
This property just follows from (13.101) . 
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Since M = 0, the expectation of M n , for n > 1, is given by E K [M n ] = E^[E^ _1 (M n )] = 
E^[M n _i] = • • • = M = 0. The mean trend of a martingale starting from zero is null. To 
estimate the fluctuations, we compute the second order moment. Setting AM,- = M J+1 — Mj 
for all j > 0, the martingale property yields E^AM,] = 0, so that, for all n > 1, 

n-1 

E K (M 2 ) =J2 e k [AM, 2 ] + 2 E * [AMjAMj] 

fc=0 0<i<j<n-l 
n-1 

= ^ E^ [AM 2 ] +2 ^ E^ \& K ( AM AM, , 

, . fc=0 0<i<j<n-l 

^' AdJ n-1 

= E k [AM 2 ] +2 E k [AME^ (AM, 

A:=0 0<i<j<n-l 
n-1 

= ^E^[AM 2 ]. 

k=0 

It remains to compute the expectation of the increments (AM^jy. 

We first prove that E^[AM 2 ] = O(hAt). Since AM fc = X k+l -~X k - E k K (X k+1 - X k ), we 
have (V# denotes the conditional variance knowing Kq, . . . , under P^) 

E^~ [AM 2 ] = V k K [X k+1 - X fe ] = E^[(X fc+1 - X k f] - (E^[X fc+1 - X k }f 
<E k K [{X k+1 -X k ) 2 ]. 

Knowing the position of the chain at time step k, the conditional probability of jumping 
is bounded by (max(a^ fc K -, 0) + max(a KktK +, 0))At/ Axk„- When jumping, the distance 

between X k and X k+ i is always bounded by Axx k - Hence, E k K [(X k+ i — X k ) 2 ] = O(hAt). 
Taking the expectation, we deduce that E^[AM|] = O(hAt). 
By ( 13. 131) . we deduce 

E/<(M 2 ) = 0(nhAt). 
By (13. 7p and (13 . 12[) and by the Cauchy-Schwarz inequality, we deduce 

Proposition 3.4. Under the assumptions introduced in the beginning of Section [3J there 
exists a constant C > 0, such that at any time n > 0, 

sup \ \u n K -u{nAt, OIIloouo < C((nhAt) 1/2 + nhAt + nAt 2 + h) exp(«nAt). 
KeT 

3.4. Interpretation by the Central Limit Theorem. This section only concerns the 
special case with constant velocity on a uniform mesh. It provides a finer result in this 
simplified case, by the use of the central limit theorem. This analysis will not be performed 
in higher dimension. We again assume that a(x) = a > and Axk is constant. We also 
reinforce the CFL condition, asking aAt < h. (This is not a restriction. When a At = h, 
the term of order 1/2 vanishes in Proposition 13.31 and the error is of order 1: this case is 
trivial.) As explained above, the order of the numerical scheme is given by the order of 
the fluctuations of the random characteristic around its mean trend. In the specific setting 
where both a and Ax are constant, the random characteristic corresponds to a random walk 
with IID increments: by the elementary theory of stochastic processes, we know that the 
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fluctuations of the walk around its mean trend are governed by the Central Limit Theorem 
(CLT in short). (See [TBI Chapter III, §3] for the standard version of the CLT and [TO], 
Chapter 2, Theorem 4.20] for the functional version in the case of a simple random walk.) 
We deduce that the fluctuations are of diffusive type, that is they correspond, asymptotically, 
to the fluctuations of a Brownian motion around the origin. This means that the random 
characteristics can be seen, asymptotically, as the paths of a Brownian motion, with a non- 
standard variance. (That is, the variance of the Brownian motion at time t isn't equal to 
t, but to a constant times t, the constant being independent of t. In our framework, the 
constant is proportional to h.) From an analytical point of view, we are saying that the 
numerical solution is very close to the solution of a second order parabolic equation: this is 
nothing but the numerical diffusive effect. 

We specify this idea. The variables (X n+1 — X n + aAt) n > are IID, with zero as mean and 
aAt(h — aAt) as variance. By the CLT, we know that 

n-1 

(3.15) (naAt{h - aAt))~ 1/2 ^(X fc+1 - X k + aAt) => W(0, 1) as n — > +oo, 

fc=0 

on each (tt, A, Pjf), K G T. The notation =^> stands for the convergence in distribution. (In 
short, for a family of random variables (W n ) n6 {o,...,+ O o} on -4.> P)> we sa Y that W n 
if Ek0(W n )] — > E[y9(H / 00 )] for any bounded continuous function (p.) The notation Af(0, 1) 
stands for the reduced centered Gaussian law. 
Plugging (13.151) into (13. 6p . we deduce that 

u n K « E K [u°(X - an At + (naAt(h - aAt)) 1/2 W)] + 0(h) 

for n large. Above, W denotes a reduced centered Gaussian random variable. The symbol 
w means that both sides are close. This point will be specified below. 
Expliciting the density of the Gaussian law, we can write 

E K [u°(X - anAt + (naAt(h - aAf)) 1/2 l^)] 

(3-16) r _ (y-anAtf dy 

J R U [XK ' K+ V) 6XP [ 2a{h - aAt)nAt^ [2va{h - aAt)nAt} 1 / 2 ' 

where xk,k+ stands for the unique point in the intersection of K and K + , i.e. the right 
boundary of the cell K when a > 0. According to [Bl Chapter 1] (with a = 0, the general- 
ization to a being trivial), we recognize the value at time nAt and at point xk,k+ of the 
solution v to the Cauchy problem: 

d(h — aAt) n0 m 
(3.17) d t v + ad x v - — — -d\ x v = 0, t > 0, i£R, 

with u° as initial condition. Finally, we can say that u\ is close to v(nAt, ■) in the cell K. 

Of course, we have to say what "close" means! This question is related to the rapidity of 
convergence in the CLT, that is the rapidity of convergence in (13.151) . The main result in 
this direction is the Berry-Esseen Theorem. (See [T5| Chapter III, §11].) In what follows, 
we use a refined version of it. (See [TBI Chapter V, §4, Theorem 14].) Denoting by F n the 
cumulative distribution function 

n-1 

\/z e R, F n (z) = F K {(naAt{h - aAt))~ 1/2 ^(X fc+1 - X k + aAt) < z}, 

k=0 
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and by $ the cumulative distribution function of the jV(0, 1) law, we have, for all z G H 

\F n {z) -${z)\ < Cn- 1 / 2 (aAt(/i-aAt)) _3/2 E^[|X 1 -X + aAt| 3 ](l + \z\)~*, 
for some universal constant C > 0. The moment of order three is given by 

E K HXx - X + aAt| 3 ] =^(h- aAt) 3 + (l - ^) (aAt) 3 
L J h h ' 

= aAt (l - ^) [(/i - aAt) 2 + (aAt) 2 ]. 

Hence, for all zeK, 

(3.18) \F n {z) - ${z)\ < C(naAt(h - aAt))' 1 ' 2 ^ 1 [(h - aAt) 2 + {aAt) 2 ] (1 + |^|)- 3 . 



By (ED, 



u n K = I u°(x KtK + - anAt + (naAt(h - aAt)) l,2 y)dF n {y) + 0(h), 
Jr 



where the integral in the right-hand side is a Lebesgue-Stieltjes integral. (See [HJ Chapter 
II, §6].) 

Assume for a while that the support of u° is compact. Performing an integration by parts 
(see (HI Chapter II, §6, Theorem 11]), we obtain (since u° is Lipschitz continuous) 



u n K = ('i i a 



At{h - aAt)) 1/2 / —(x KiK + -anAt + {naAt{h-aAt)) 1/2 y)F n (y)dy + 0(h), 
Jr dy 

Plugging (I3.18P in this equality, we obtain 

f du° 

\u n K - (naAt(h - aAt)) 1/2 / — (x K>K + - anAt + (naAt(h - aAt)) 1/2 y)$(y)dy\ 

Jr dy 

< Cnh- 1 [(h - aAt) 2 + (aAt) 2 ] f (1 + \y\)~ 3 dy + 0(h). 

Jr 

Performing a new integration by parts and then a change of variable, we see that the left- 
hand side is equal to \u\ — v(nAt,x KyK +)\, where v is the solution of the Cauchy problem 
f)3.17p . Using a standard troncature argument, we can easily get rid of the assumption made 
on the support of u°. We finally claim 

Proposition 3.5. Assume that a(x) and Axk are constant and that u° is n-Lipschitz con- 
tinuous. Then, there exists a constant C > such that, at time any time n > 0, 

sup \ \uk - v(nAt, OIIl^/o < Ch, 

where v stands for the solution of the Cauchy problem (I3.17P with u° as initial condition. 

It is remarkable that the bound doesn't depend on (n, At). The result is still true for 
aAt = h. (See Proposition 13.31 ) 
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4. Principle of the Analysis in higher dimension. Application to a Simple 

Case 

We here present the basic ingredients for the analysis in dimension d greater than two. 

As above, the velocity field a is assumed to be bounded and ft-Lipschitz continuous. The 
regularity of the initial condition u° will be specified below. 

The characteristics of the transport equation are still denoted by (Z(t,x)) t >o, x G IR d , see 
Equation (TO) . 

As in [13], the cells are assumed to be uniformly non flat, i.e. they satisfy, in a strong 
sense, the converse of the isoperimetric inequality: 

(4.1) 3a > 0, VK G T, ^ \K n L\ < ot\K\hr x . 

By the standard isoperimetric inequality, this is equivalent to the existence of (3 > 1 such 
that \K\ > (3~ l h d and J2l~k \K D L\ < ph^ 1 for all K G T. 
The mean velocity on a cell K is denoted by 

ax = {K^ 1 / a(x)dx. 
Jk 

The constants "C" and "c" below may depend on ||a||oo ; « ; P, k> and d. As in dimension 
one, they are always independent of At, of h, of the current time index n and of the random 
outcome u. 

4.1. Stochastic Representation of the Scheme. The expression of u 1 ]^ 1 given by the 
upwind Scheme is (see (12.21) ) 

n+i (a K ,L,n K ,L}At\K HL\ ( ^ (a K>L , n KtL )At\K n L\ \ n 

u K =- 2_. \k\ ^ + 1 + 2^ \k\ u ^ 

where K~ = {L ~ K, (cbK,Li ^k,l) < 0}, u K being given by 

u° K = IK]- 1 [ u°(x)dx. 
'k 



In this framework, the CFL condition has the form 

\/K e r - ( aK > L > UK > L } At \ KnL \ < i 



Lei<- 



and is assumed to be satisfied in all the following. 
As in dimension one, the coefficients 



(a KjL ,n K>L )At\K n L\ 
Pk,l = VrF\ tor L e K 



Pk,k = 1 + 



(aK,L,n K ,L)At\K n L\ 
\K\ 



LeK- 

p K)L = for L G T \ (K~ U K) , 

can be interpreted as the probability transitions of a Markov chain with values in the set 
of cells. Again, we can find a measurable space (ft, A), a sequence (K n ) n > Q of measurable 
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mappings from (Q, A) into T as well as a family (Pe-)xgt of probability measures on (Q, A), 
indexed by the cells, such that, for every cell K e T, (K n ) n >o is a Markov chain with K as 
initial condition and {j>k,l)k,l&t as transition probabilities. As in dimension one, the chain 
(K n ) n > Q goes against the velocity field a: for n > 0, either X n+ i is equal to K n or i^ n+ i 
belongs to K~. Similarly, for n > 1, either K n _i is equal to K n or K n _i belongs to K+. 

Following the analysis performed in dimension 1, we can prove the backward Kolmogorov 
formula: 

Theorem 4.1. Under the above notations, the numerical solution u K at time n and in the 
cell K has the form: 

< = E X [<J. 

Due to the strong similarity with the one-dimensional frame, we do not repeat here the 
proof. 

In what follows, for K e T, we denote by E^- the conditional expectation Ex[-|-K"o, • • • , K n ]. 

4.2. Random Characteristics. To follow the one-dimensional strategy, we have to as- 
sociate a sequence (X„) n > of (random) points with each path of the Markov chain. In 
dimension one, the point X n is defined as the entering point in the cell K n . Two points may 
play this role in the higher dimensional setting: either the barycenter of the entering face 
K n _i n K n or the barycenter, with suitable weights, of all the possible entering faces in the 
cell K n . For a given cell K, we thus define xk,l as the barycenter of the face K n L if K 
and L are adjacent: 



cc doc ■ 

KnL 



dK,j = for J e K\ 



(4.2) x k< l = \KnL\- 1 [ 

J j- 

and Ck as the following convex combination of (xk,j) jeK+' 

(4-3) e K = MT q K ,j £ 

\JeK+ J JeK+ 

with 

(a>K,j,n K ,j)At\K n J 
\K\ 

(4 4) . ( a K,j,n KtJ )At\K n J\ 
' qK ' K = 1 ~ 1^ \k\ ' 

J&K+ 1 1 

q KJ = for J e T \ {K + U K) . 

We emphasize that the {qK,.j).j&T are, at least in a formal way, the weights associated with 
the scheme for the velocity —a. We will specify this correspondence below. We also notice 
that ck might be outside K if K is not convex. This does not matter for the analysis. 
With these notations at hand, we define the random characteristic X n as 

Xq = ex , 

(4.5) X n = e Kn if K n = X„_i, n > 1 

X n = X Kn _ u K n if Kn ^ Kn-l,n > 1. 
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This definition is quite natural. When n = or K n = K n _i, n > 1, X n is chosen as a 
remarkable point of the cell K n , independently of the past before n. When K n ^ K n _i, the 
choice of X n expresses the jump from K n _i to K n . 

4.3. Green's Formula. The result given below explains why the barycenters of the faces 
are involved in our analysis. The main argument of the proof relies on the Green formula, 
which has a crucial role in the whole story, as easily guessed from the specific form of the 
transition probabilities. (See also [U, Proposition 3.1].) 

Proposition 4.2. Consider a cell K . For any point xq in the convex envelope of the cell 
(we say convex envelope because of ex, defined as a barycenter) , 

a K At } Pk,l(xk,l -xo) + ^ qK,j(x K>J - x ) + O(hAt). 

L£K~ J&K+ 

Proof. For an index 1 < i < d, the Green formula, see [HI Chapter 3, (3.54)], yields (a^ 
and Xi stand for the i th coordinate of a and x). 

(4.6) / a,i(x)dx = - (xi - (x )i)div(a)(x)dx + y~] / (x^ - (x )i){a(x), n K ,L)dx. 

JK JK Lr^K J KnL 

The left-hand side is equal to \K\ (ax)i- By the regularity of a, the first term in the right-hand 
side is bounded by 0{h\K\). Similarly, the last term in the right-hand side writes 

L~K jKnL 

= ^2(a KtL ,n KiL ) / {xi-{x )i)dx+^2 / fa - (x )i)(a(x) - a K)L ,n KjL )dx 

L~K jKnL L~K jKnL 

= ^ {aK,Li n K ,L) \ (xi - (x )i)dx + 0(h 2 )^2\K f] L\. 

L~K jKnL L~K 

With ffl~T]) and (fl~2j) at hand, we deduce that 



(4.7) L~K jKr < L 



t . K JKnL 

J2 ( a K,L,n K ,L)\K n L\ [x KiL -x ]. + 0{h\K\). 



L~K 

From (14. 6p and (14. 7p . we claim 



a, (a K ,L,n KtL )At\K nL\ f , n/ , A . 

a K At = — [x k ,l - xo) + O(hAt). 



L~K 



This completes the proof. □ 
The following corollary is the multi-dimensional counterpart of Proposition 13.21 

Corollary 4.3. For any starting cell K G T (so that we work under F^) and any n > 0, 

El [X n+1 - e Kn ] = -a(X n )At + O(hAt). 
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As a simple consequence, when the cell K n has only one possible entering face (entering 
means entering for the random characteristic), as it is the case in dimension 1 with a > 0, 
ex n = X n and the conditional expectation of the mean displacement X n+ \ — X n knowing 
the past before n is driven by —a, as stated in Proposition 13.21 

Proof. For K G T and n > 0, 

(4-8) E™ [X n+1 - e Kn ] = PK n ,L {x Kn ,L - e Kn ) . 

L&K- 

(Indeed, if K n+ % = K n , X n+ i = ex n -) Applying Proposition 14.21 with xq = 6k„, we obtain 
(4.9) a Kn At = - VK n ,L(x Kn>L - e Kn ) + ^ lK n ,j{x Kn ,j ~ e Kn ) + O(hAt). 

L€K- JeK+ 

By the very definition of e^„ , the second sum in the above right-hand side is zero. Identifying 
(14. 8 p and (14. 9p . we complete the proof. □ 

4.4. Set of Problems. Keeping the one-dimensional strategy in mind, we understand that 
the whole problem now consists in estimating the gap e^ n — X n for n > 1. (For n = 0, it is 
zero.) 

As said above, ex n — X n vanishes when the cell K n admits only one entering face. 
Unfortunately, there is no hope to obtain a similar result, or an estimate of the form 
ex n — X n = O(hAt), for a cell K n of general shape. 

The reason is purely geometric and is well-understood on a two-dimensional triangular 
mesh. Except specific cases, the cardinality of K + , K being a triangle, is either one or two. 
(See Figured] below.) 




Figure 1 . Example for the gap e^ n — X 5 



In this picture, the velocity a is constant, oriented from the top to the bottom. The 
triangles J\, Ji and L have only one entering edge and the triangle K has two entering 
edges. By Corollary I4.3[ the mean displacements in the triangles J\, J 2 and L are driven by 
—a. The story is different for K. 
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The position of ex is determined in the following way: e# belongs to the segment 
[ x k,Ji, x k,j 2 ] an d, by Proposition 14. 2\ xk,l — &k is parallel to a (the term O(hAt) in Propo- 
sition [472] is zero since a is constant). Assuming that at time n — 1, K n ^\ = J\ and that, at 
time n, K n = K, X n is the middle of the edge K PI J\. We understand that, in this case, 
ex n — X n is of the same order as h. 

The key point in our analysis follows from a simple observation: when reversing the velocity 
field a in Figure dj i-e. when changing a into —a, the triangles with one entering edge turn 
into triangles with two entering edges and the triangle with two entering edges turns into a 
triangle with one entering edge. Said in a very naive way, the bad triangle K for a becomes 
a good triangle for —a. 

Here is another simple fact: the deterministic characteristics for —a correspond to the 
characteristics for a, but reversed in time. 

These two observations lead to the following idea: in what follows, we hope to get rid of 
the remaining gaps (e^ n — X n ) n >i by reversing the random characteristics. 



4.5. Reversing the Markov Chain in the Divergence-Free Setting. As announced 
above, the weights (qK,j)K,jeT in (14.41) correspond, at least in a formal sense, to the scheme 
associated with the velocity field —a. We say "formally" because the CFL condition for the 
field —a may fail: as shown below, the term qx,K m ay be negative. 

However, when a is divergence- free, the CFL condition holds for —a and q is the analogue 
of p for the field —a: 

Proposition 4.4. Assume that a is divergence- free. Then, 

so that qx,K = Pk,k for every cell K . In particular, the CFL condition holds for —a and q 
is a Markovian kernel, i.e. qx,j > for all K, J G T and YljerQK'J = 1- 

Proof. The proof follows again from Green's formula. Indeed, for every cell K, 

(4.10) / dxv(a)(x)dx = 0=y2\KnL\(a K ^,n KtL ). 

Jk l~k 

The terms in the right-hand side are equal to — pk s l\K\/ At if L G K~ and qx,L\K\/ At if 
L G K + . This completes the proof. □ 

In the rest of this section, we assume a to be divergence-free. The point is to understand the 
connection between the original Markov chain (i^ n ) n >o, associated with a, and the Markovian 
kernel {qK,j)x,j&T- In what follows, we show that, under suitable conditions, the random 
process obtained by reversing the chain (K n ) n >o is a Markov chain with (qK,j)K,jer as 
probability transitions. 

Reversing a Markov chain is a standard procedure in probability theory. (See e.g. [T5J 
Section 1.9].) The reversed process is always a Markov chain, but the transition probabilities 
may be highly non trivial and do depend, in almost every case, on time. Anyhow, the law 
of the reversed chain is easily computable when the chain is initialized with an invariant 
probability. 
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Assume for the moment that the chain (i^ n ) n >o admits an invariant probability u, i.e. 
there exists a probability u on the set of cells such that 

VKeT, v(K) = J2"(J)pj,k, 

JeT 

and pick the initial cell of the Markov chain randomly, with respect to the probability measure 
v. We then denote the corresponding probability measure on (fl, A) by F u as we did for a 
deterministic starting cell. Under P^, Kq is random and its distribution is given by v itself. 
The measure P^ may be decomposed along the measures (F^KeT according to the formula 
F u = ^2 Ker v(K)Fk. By invariance of z/, K n follows, for each n > 0, the law z/, i.e. 

Vn > 0, VK e T, F v {K n = K} = v(K). 

(See [151 Section 1.7].) Denoting by (Kq , . . . , K N ) = (Kjy, . . . , K ) the reversed chain from 
time iV to time 0, we have, for K,JeT and < n < N — 1, 

(4.11J F v {K n+x = J \K n = K) = — — — = ——p JK: 

F u {K N _ n = K\ v(K) 

so that the reversed chain is homogeneous (i.e. the transition probabilities don't depend on 
time). 

Since the velocity field a is divergence-free, the transport equation conserves the mass with 
respect to the Lebesgue measure. In this framework, the Lebesgue measure is invariant for 
the Markov chain (i^ n ) n >o: 

Proposition 4.5. Assume that a is divergence- free. Then, the Lebesgue measure is invariant 
for the Markov chain, i.e. 

VKeT, \K\ = J2\J\PJ,k- 

JeT 

Proof. For a given cell K, we have 



^2\J\pj,k= ^2 { a J,K,n JjK )kt\J (1K\ + \K\ + ^2(a JiKi n JtK )At\JnK\. 

By flEDD, 



(4.12) 

Jer JeK+ j&k 



^2 \ K n J \ ( a J,K, n J>K ) = 0, 



J~K 

since a is divergence-free. Returning to f)4.12p . we obtain J2jer \-Avj,k = \K\. □ 

4.6. Analysis in the Divergence-Free and Periodic Setting. The Lebesgue measure 
is not of finite mass on Mr, but it is a probability measure on the torus R d /Z d . 

To explain how we make use of time reversal, we thus assume, for the moment, that the 
problem is periodic, of period one in each direction. (Since the transport is of finite speed, 
this is not a big deal. Anyhow, the proofs of the main results of the paper, given in the next 
section, are performed without any periodicity assumption. The current paragraph is purely 
pedagogical.) This means that both the velocity a and the mesh are periodic, of period one 
in each direction of the space. 

As a consequence, we can see the Markov chain (K n ) n >o as a Markov chain with values in 
the set T/Z d , i.e. in the space of classes of cells for the equivalence relation induced by the 
periodicity. The probability of jumping from the class of the cell K to the class of the cell 
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L is given by the rate Pk,l- by periodicity, this rate doesn't depend on the choices of K and 
L. 



Proposition 4.6. Assume that a is divergence- free and that both a and the mesh are periodic 
of period one in each direction of the space, then the Lebesgue measure on the torus induces 
an invariant probability for the Markov chain, i.e. 

VKeT/z d , \K\ = \ j \pj>k- 

Jer/z d , J~K 

(For J G T/Z, d , we denote by \J\ the common volume of all the cells of J.) 
In particular, by (14.111) . the reversed chain has the following transition probabilities 

P M {iC +1 = J\K N n =K}= q K>J . 

(In the above equality, /i denotes the Lebesgue measure on the torus. Under F^, the starting 
class of cells of the chain is chosen according to the Lebesgue measure.) 

Proof. To check the last assertion, we note that, for two adjacent cells J and K, p JK > 
if and only if K G J", that is J G K + . Hence, pj^k > -v=> qx,j > and in this case 

( 4 - 13 ) <1k,j = T^\PJ,k, 

so that the relation is always true (even if one term vanishes). □ 
We recall the interpretation of Proposition 14.61 the reversed chain (seen as a chain with 
values in T/Z d ), when the chain (K n ) n > is initialized with the Lebesgue measure, is nothing 
but the chain associated with —a. 

Before applying Proposition to the analysis of the numerical scheme, we have to specify 
the construction of the random characteristic (X n ) n > in the periodic setting. In this case, 
(X n ) n >Q is seen as a path in M. d and not in M d /Z d . This amounts to consider a sequence 
of representatives (K n ) n > for the chain (K n ) n > . We choose K as the only representative 
Kq of K such that e K o has coordinates in [0,1), that is K = Kq (and X = e K o). For 

any n > 0, K n +i is the unique representative of K n+ i such that K n +i G {K n )~ . With this 
sequence of representatives at hand, we can build up the sequence (X n ) n > according to 
(14. 5p . We are now in position to estimate the gap e^ n — X n , n > 1: 

Proposition 4.7. Under the assumptions of Proposition pJT6| for any N > and any n G 
{0,...,N}, 

^[e kn -X n \K n ,...,K N ] =0, 

E M denoting the expectation under P M . In particular, under F^, the sequence (M^)o< n <jV; 
given by 

_ _ JV-l 

K = Mf = 0, M N n = ( e K k ~X k ), 2<n<N, 

k=N-n+l 

is a martingale for the backward filtration (a(K]y_ n , . . . , K^) = o-(Kq , . . . , -ft',^))o<n<A r - 

Proof. We can assume that 1 < n < N, since e ko — X = 0. We then emphasize that K n 
isn't measurable with respect to a(K n , . . . ,K^). This is the main difficulty of the proof. 
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Indeed, K n depends both on K n and on the initial representative of K . Anyhow, the 
difference eg — X n doesn't depend on the representatives chosen for K n -\ and K n . 

Indeed, recalling that K® is the unique representative of K n such that e#o has coordinates 
in [0, 1) , we can always write 

( 4 - 14 ) e K n ~ X n= ( e K n - X K n ,j) 1 {k n . 1 =J}= ( e K° - X K O tJ )l {Kn _ lBj} . 

Je(k n )+ J£(K°)+ 
(Above, K n _i 3 J means that J is a representative of K n _i.) Hence, 

E v[ e k n - X n \K n , . . . , K N ] = Y (eK»-x K o :J )F^{K n ^ 1 3J\K n ,...,K N }. 

Je(K°)+ 

The probability F^K^ 3 J\K n ,...,K N } is also F^{K N N _ n+1 3 J\K%, . . . , K N N _J. By 
Proposition 14.61 it is equal to q K o j. By (|4.3p . we deduce that E^e^ — X n \K n , . . . , K N ] = 0. 
By ( 14.141) . we know that efc n 1 — A n+1 is o~(K n , . . . , K N ) -measurable for < n < N — 1. We 
easily deduce the martingale property. □ 
Following (I3.12p . (13.131) . (I3.14p . we complete the analysis by estimating, under P M , the 
fluctuations of the random characteristic around the deterministic characteristic. 

Proposition 4.8. There exists a constant C > such that, for any N > 1, 

^[\X N - Z(NAt,X )\] < C[(NhAt) 1/2 + NhAt + iVAt 2 ] exp(nNAt). 

Proof. Following the proof of (13.121) and keeping the identity e^ o — X = in mind, it is 
sufficient to focus on 

jV-1 N-l N-l 

X N -X + J2 a{X k )At = ^(A fc+1 - e kk + a(X k )At) + ^ e k k ~ X ^)- 

k=0 k=0 k=l 

It is plain to see that Corollary 14. 31 is still true under P^ (and not under F K ), so that 

N-l N-l N-l 

X N -X + J2 <X k )At = i x *+i ~ tk k - Ej(A fc+1 - e kk )) + ~ **) + O(NhAt). 

k=0 k=0 k=l 

(Above, Ej[-] = E M [-|K , . . .,K k \.) Setting, for all < n < N, M n = Yllli^k+i - e kk - 

E* 
n 

write 

N-l 



Fr(X k+ i — X k )) (with M = 0), and using the notation introduced in Proposition 14.71 we 



X N - X + J2 a ( x k)^t = M N + M N N + O(NhAt). 

k=0 

We let the reader check that (M n )o< n <7v is a martingale with respect to the (forward) filtra- 
tion (a(Ko, . . . , K n )) < n <N. By Proposition 14.71 (M^) < n <7v is a martingale with respect to 
the backward filtration. 
Following (13.131) . we have 

N-l _ N-l _ 

M\ M N\ 2 ] = J2M\ M n + i-M n \ 2 }, e,[\m n n \ 2 ] = Y E »[\K +1 -K\ 2 ]. 

n=0 n=0 
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Following (13.141) . we have, for all < n < iV — 1, 

E4|M n+1 - M n \ 2 ] < E^[\X n+1 - e k f] 

, x <M PKu,L\ x K n ,L-^K n \ 2 } 

(4-15) Le( * n) - 

< h 2 sup Y Pk,l < llallooSupfliTl" 1 ^ \KnL\]h 2 At. 

K L~K K L~K 

By (14. II) . we deduce E^jjMn+i — M„| 2 ] < aHaHoo/iAt. By a similar argument, we obtain 
E M [|M^ +1 - M N n \ 2 } < aWaWovhte. We deduce 

E^[\M N \ 2 ] =0(hNAt), E^[\M N N \ 2 ] =0(hNAt). 
We then complete the proof as in the one-dimensional case. □ 

4.7. L 1 -Error in the Divergence- Free and Periodic Setting with u° Lipschitz con- 
tinuous. As a by-product, we obtain the following estimate for the error of the numerical 
scheme when u° is Lipschitz continuous: 

Theorem 4.9. Assume that the hypotheses of the beginning of Section^ are satisfied. As- 
sume moreover that a is divergence-free and that both the velocity a and the mesh are periodic 
of period one in each direction of the space. Assume also that u° is n-Lipschitz continuous. 
Then, there exists a constant C > such that 



E 



' \ux - u(NAt,x)\dx < C((hNAt) 1/2 + hNAt + NAt 2 ) exp(/dVAf). 

K 



Proof. Consider a cell K. We know that is nothing but 



u 



E k [u° Kn ] = E K [u°(X N )]+0(h). 



Moreover, for all x £ K, \Xq — x\ < h with probability one under P#. By stability of the 
solutions of (Oj) . \Z{NAt,x) - Z{NAt,X )\ < hexp{nNAt) under F K . We deduce 

Vx £ K, u{NAt,x) = u°(Z(NAt,x)) = E K [u° (Z(NAt, X ))] + O (h exp(nN At)) . 

Hence, 

J \u% -u{NAt,x)\dx < k Y \K\E K [\X N -Z(NAt,X )\]+O(hexp(KNAt)) 

KeT/Z d K Ker/z d 

< kE^ [\X n - Z(NAt, X ) \]+0(h exp(nNAt)) . 
This completes the proof. □ 

5. Analysis in the General Setting 

We now turn to the general case and analyze the error of the numerical scheme, both in 
the L 1 sense and in the L°° sense, u° being respectively of bounded variation and Lipschitz 
continuous. We thus forget the periodic setting and the divergence-free condition. Anyhow, 
for technical reasons explained below, the time step At is required to be small when the 
divergence of a is large, i.e. 

(5.1) 3t7£(0,1), ||div(a)|| 00 At<l-r / . 
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We keep the notations of Section HI To simplify the form of the final bounds, we assume 
that h < 1 and that At < 9h for some 9 > 0. The constants C, C and c below may depend 
on the parameters specified in Section HI on rj and on 9. The values of these "constants" 
may vary from line to line. 

5.1. Strategy. In dimension one, we were able to analyze the error in the L°° sense by 
investigating the distance between X N and X — At J2n=o a (X n ) for any arbitrary initial cell 
K of the mesh. In the previous section, the result was given in the L 1 norm since the chain 
was initialized with the Lebesgue measure. In what follows, the approach is halfway. 

The idea is the following. We pick up the starting cell of the chain with respect to the 
Lebesgue measure among the cells included in a ball of radius h 1 ! 2 and centered at the origin 
(or at any other arbitrarily prescribed point). The bounds we then obtain for the error of the 
numerical scheme hold in the L 1 sense, but locally in the ball. In other words, we manage 
to bound 

/T d/2 sup f \u N (y) -u{NAt,y)\dy 

xm d JB(x,h 1/2 ) 

by a constant times h 1 / 2 (up to remaining terms in N and At) , u° being Lipschitz continuous. 
(In the above expression, u N (y) stands for u 1 ^ when y belongs to the interior of K.) By an 
approximation procedure, we deduce that the scheme is of order 1/2 for the (global) L 1 norm 
when u° is of bounded variation. 

Actually, we can perform the same analysis by replacing the local L 1 norm by the local LP 
norm, p being greater than one. We then derive that, for any small positive e, the scheme is 
of order 1/2 — e for the L°° norm when u° is Lipschitz continuous. Moreover, by translation, 
it is sufficient to prove this estimate for x — 0. 

For all these reasons, the quantities of interest are 

N-l 

Qp = h d ' 2 [\X N -X + AtY, a(X n )\ p ] , N,p>l, 

KeT n=0 

where % stands for the set of cells K such that |e#| < h 1 ! 2 . By ( 14. ip . the cardinality of % 
is bounded by Ch~ d l 2 for some positive constant C . 

For given N, p > 1, we are going to decompose Qp along all the possible paths of the 
chain. To do so, we distinguish the paths according to their fluctuations around the velocity 
field —a. We write 

n-l 

Qp < h id+p)/2 J2( k + ^K{kh 1 ' 2 < sup \X n -X + AtJ2^)\ <(k + l)/i 1/2 }. 

js^n- i Kn<N . n 

KgTq k>0 — — t=0 

To estimate the above right-hand side, we will use a time reversal argument, as in Section 
HJ Therefore, we are mainly interested in the location of the arrival cell on the event 
{kh 1 ' 2 < su Pl < n <jv \X n - X + AtJ2^a(Xi)\ < (k + l)^/ 2 }. Thus, for any k > 0, we 
denote by T£ the set of cells such that there exists an iV-tuple ( J , . . . , Jjv-i), with 
Jo G %, satisfying, for all i G {0, . . . , N — 1}, either J i+1 = Ji or J i+1 £ Jf, and 

n-l 

(5.2) kh^ 2 < sup \y Jn _ uJn -ej +AtJ2<yj i -uJ i )\<(k + l)h 1/2 , 

l<n<N ^ 

23 



with y JuJl+1 = ejlf J { = J i+1 and y Ju j i+1 = x JiiJi+1 if J i+1 G J { (and yj_ 1)Jo = e Jo ). 
Under Pjf, for K G %, we have -KOv £ 7^ on the event {kh 1 ^ 2 < sup 1<n<Ar \X n — X + 

Atz:: 1 *(x*)\<(k + i)h 1/2 }. 

Thus, 
(5-3) 

xP K {tf/ 2 < sup |X„-X + At Va(X;)| < (A; + l)/i 1/2 ,/i~jv = L}]. 

l<n<iV ^0 

In the next lemma, we estimate the cardinality of 7^. Because of the translation action 
of the transport equation, $[7^] is of the same order as |}[7o]: 

Lemma 5.1. There exists a constant C > suc/i t/jat ? /or all k > 0, 

%[T k N ) < Cex V {CNAt){k + l) d h- d/2 . 

Proof. We fix k > and we consider a sequence of cells (Jo, • • • , Jn), Jo £ To, satisfying 
O. Setting . . . ,2/^) = (e Joi yj 0iJl , . . .,yj N _ u j N ), 

n-1 



sup \y n -y° + AtJ2 a (v l )\ < {k + l)^' 2 . 



Kn<N 



i=0 



Plugging the characteristic of the transport equation (see (11.21) ). we deduce (recall that h < 1 
and At < 0/i by assumption) 



n-1 



sup \y n - Z(nAt, y°) + At V [a(j/*) - a(Z(iAt, i/ ))] 

<n<N ^ 

<{k + l)h 1/2 + CNAt 2 <{k + l)h 1/2 {l + CNAt), 



for some constant C > 0. By the Lipschitz property of a and the Gronwall lemma, it is plain 
to deduce (up to a new value of C) 

\y N - Z{NAt,y°)\ < C exp{CNAt)(k + l)h 1/2 . 

Since \y°\ < h 1 / 2 (J G To), we deduce, by stability of the solutions to (11.21) . that 

\y N - Z{NAt,0)\ < Cexp{CNAt){k + l)h 1/2 . 

Every point x G Jn satisfies the same property since the diameter of Jn is bounded by h. We 
deduce that there exists a constant C such that Jn is included in the ball of center Z(NAt, 0) 
and of radius Cexp(CNAt)(k + l)h}^ 2 . In other words, all the cells in are included in 
this ball. Up to a modification of C, the volume of the ball is Cexp(CA r At)(A; + \) d h d l 2 . 
Since the volume of a given cell is greater than (3~ l h d (see Assumption (14. ip ). the cardinality 
of T£ is bounded by C exp(CNAt)(k + \) d hr d l 2 for a new value of the constant C. □ 
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5.2. Application of Section [4], In light of Section HJ we introduce the following decompo- 
sition 

N-l 

Xn — Xq + At a(X n ) = Sn + Rn, 

n=0 

with 

n— 1 n— 1 

S = R = 0, S n = Y,l X i+i ~ e Ki + &ta(Xi)} } R n = ^[e Ki -Xj, n > 1. 



i=0 i=0 



On the event {fc/i 1/2 < swp x < n < N \X n - X + At^vTc i a PQI < ( fc + l ) hl,2 }i we have 
su Pi<n<Af I'S'nl > {k/2)h 1 / 2 or sup 1<n<7V |i? n | > (k/2)h 1 / 2 . By (15.31) . we obtain 

rWQf<rV V (A; + l)"P*-{ snp |S n |>V 2 ,i^ = L} 

(5.4) 

4- V V V (/,■+! )'■¥>,, i sni) /?.. > y 



+ EE E( fc+1 ) P M su p i^ni>|^ 1/2 ,^ = ^}. 



fc>0 KeTo Le7 N °< n < N 



We wish to apply Corollary 14.31 to treat the first term in the above right-hand side. Since 
11 [73] < Ch- d / 2 , we have 

EE E( fc + 1 ) P M su p \s n \>h i / 2 ,K N = L} 

T^n t77Z- — v 0<n<iV 2 

k>0 K&T LeT^ ~ ~ 

(5.5) <Ch- d / 2 J2(k + iy supF K { sup \S n \>^0 2 } 



< ch- d ' 2 



h ^NAt(i + cNAty + + i) p bM-cm;) + exp (-c^)] 

fc>0 



the last line following from 

Lemma 5.2. There exists a constant C > such that, for k > CNAt h 1 / 2 , 

Proof of Lemma l5.2l We fix the starting cell K e T. (We thus work under Pa-.) Following 
the proof of Proposition I4.8[ we introduce the sequence 

n-l 

M = 0, M n = - - E^[X i+1 - e*-J), n > 1. 

i=0 

It is a martingale with respect to the filtration (a(K , . . . , i^ n ))„>Q (under P^r). By Corollary 
14.31 we have 

sup \M n - S n \ < CNAth, 

0<n<N 

for a positive constant C. Thus, on the event {sup 0<n<Ar l-S^I > (k/2)h 1 ^ 2 }, sup 0<n<N \M n \ > 
(k/A)^ 2 or CNhAt > (k/^h 1 / 2 . The latter is impossible if k > ACNh 1 / 2 At~ We deduce 
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that 

sup \S n \ > k -h xj2 } <¥ K { sup (iky^-Zi 1 / 2 }, 

0<n<N 4 0<n<N 4 

for k > ACNAth 1 ^ 2 . As in (14.15p . the conditional variance of the martingale (M n ) n >o may 
be bounded by CAth for a possibly new value of C: 

Vn > 0, E™ [|M n+ i - M n \ 2 ] < CAt h. 

Moreover, the jumps of the martingale are bounded by h, i.e. |M n+1 — M n \ < h for all n > 0. 
Applying Proposition 16.11 given in Annex to (/i _1 M n ) n >o, we obtain 

F K { sup \M n \ > k -h 1 ' 2 } < C[exp(-^|^) +exp(-^_)]. 



0<n<7V 



(Pay attention, the value of v in Proposition 16.11 is v = CAth 1 since we divide (M n ) n > by 
h.) 

5.3. Time Reversal. To treat the gap term i n (15.4|) . we are to reverse the chain (K n ) n > 
as done in Section H] and then to compare the law of the reversed chain with the law of the 
chain associated with —a. 

Because of (I4.10p . we emphasize that the CFL condition may fail for —a, so that we cannot 
associate a Markov chain with the weights {qK,j)K,j&T- The following proposition says how 
to modify them to obtain a Markovian kernel: 

Proposition 5.3. Set 

\/K e T, 5 K = IK]- 1 [ div(a)(x)dx. 

Then, under Condition (15. ip . the kernel 

lK)J = (1 + 8 K At)- l q K ,j for J G K + , 
1k,k = 1 - ^ (1 + 5 K Aty 1 q KJ , 

J&K+ 

7XiJ = 0for J eT\(K + UK) , 

satisfies "/k,k — (1 + &k^)~" x Vk,k- In particular, it is Markovian, i.e. jk,j > for all 
K,J G T and Y.jer "1k,j = 1- 

By (15. ip . we emphasize that jk,j < V~ 1( 1k,j for J G and 7a-,x < Tj~ l pK,K- 

Proof. Consider a cell K. By (15. ip . we have 1 + 5 k At > 0. We compute "/k,k- By 
integrating by parts the expression of 5k, as in (14.101) . we get 

- ^2 Pk < l + qK ' J = S K^t, 

that is Pk,k~ l+(l+5xAi)(l— Jk,k) = We deduce that ^k,k = (l+^Ai) -1 ^,^- D 

The chain associated with the kernel 7 is denoted by the pair ((Qk)k€T, (rn)n>o) : (T n ) n > 
is a sequence of measurable mappings from (Q,A) into the set of cells and (Qk)kgt is a 
family of probability measures on (Q,A), such that (T n ) n >o is, under Q^, a Markov chain 
with X as initial condition and 7 as kernel. The expectation under Q^- is denoted by E® . 
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The link between the chain (T n ) n > and the reversed chain (K n = KN-n)a<n<N is given 

by 

Proposition 5.4. There exists a constant C > 1, such that, for any pair of cells (K, L) and 
any function \& : T N+1 — ► M +; 

C- 1 [1 - Hdiv^lUA^Eg^^o, • • • , T n )1 {Tn=k} ] 

<E k [^(K^,...,K n )1 {Kn=l} ] 

< C[l + Hdiv^llooAt^Eg^ro, • • .,T N )l {rN=K} ]. 

Of course, Proposition 15.41 is weaker than Proposition 14.61 Above, we are just able to 
compare the law of (T , . . . ,Tn) under the measure 1{t n =k}-Ql with the law of (Kn, . . . , K Q ) 
under the measure 1{k n =l} • IV- They are equivalent and the resulting density is bounded 
from above and from below by positive deterministic constants. 

Proof. Without loss of generality, we can assume that \I/ = l(j 0) ...,j JV ), with (J , ...,J N ) e 
T N+1 , J = L and J N = K. Then, 

E K . . . , K n n )1 {Kn=l} ] = F Jn {K% = J , . . . , K* N _ X = J N -i, K n n = J N } 

( 5 - 6 ) = W Jn {K = J N , K l = J N _ 1 ,...,K N = J } 

PJn,Jn-iPJn-i,Jn-2 ■ ■ m PJi,Jo m 

By Proposition 15.31 and (14.131) . we know that 

7j„,j„+i = i t wTTJT A+\ PJn+uJ n > < n < N - 1, 
\J n \{l + o Jn l\t) 

even if J n = J n+ \. Plugging this relationship in (15. 6p . we deduce that 

- JV-1 17 1 

E k [*(K»,...,K n )1 {Kn=l} ] = [n( 1 + ^» A *)]|TT7Jb,Ji---7J w - 1 ,J w 



n=0 
N-l 



[l[(l + 8 Jn At)] J^E? [*(T , • • • , r N )i { r N =K } ] • 



i - AT 

n=0 

-lud ^ I Tl ^ ud 



By (O) and (EEJ (recall that (O) implies /T 1 /^ < |J| < h d for all J G T and for some 
(3 > 1), we complete the proof. □ 

5.4. Analysis of the Gap. Using the previous subsection, we analyze the term Rn in (15. 4p . 
For N > 1, we emphasize that 

n-l 

sup \R n \= sup |^[(e^ -x^^Jl^^}]! 



Kn<N Kn<N 



A' 



<2 sup |y)[(ejf. -arKi-LJcJl^.!^}]! 

Kn<N — 

— — i=n 

n-l 

= 2 sup E^e^-x^^ )1^ ^ ]|. 

Kn<N n K i K i ' K i + 1 i"i+l^ A i 1 

— — J=0 

27 



Keeping Proposition 15.41 in mind, we define the analogue, but for the chain (r n ) n > , that is 

n-1 

H = 0, E n = ^[(e r , -xr,,r l+1 )l{r l+1 ^rj], n > 1, 

8=0 

By Proposition 15.41 we deduce that, for any starting cell K, for any terminal cell L and for 

any k > 0, 

(5.7) P A -{ sup \R n \>%h 1 ' 2 ,K N = L}<C(l + At) N Q L { sup \~ n \>\h l l\T N = K}. 

l<n<N * l<n<N 4 

Here is the main argument of the analysis. 

Proposition 5.5. For a given cell L e T ', the process (S n ) n > is a martingale under Ql 
with respect to the filtration (cr(r , . . . , r„)) n > . 

Proof. The proof is quite obvious. For each n > 0, S n is measurable with respect to 
cr(r , • • • ,T n ) and 

E ?[ e r„ - a;r n+1 ,rjr , . . . ,T n ] = ^ 7r n ,j(er„ - xj,r n ) 

Jer+ 

= (1 + 6 Tn At)' 1 £ gr„,j(e r „ - x J>Tn ) = 0. □ 
Jer+ 

Following the proof of Lemma 15.21 and making use of fl5.ll) to bound the conditional 
variances of the increments of (S„) n > , we deduce 

Lemma 5.6. There exists a constant C > such that for any k > 

supQ L {jupJS n | > h 1 ' 2 } < C[exp(-^^) +exp(-^_)]. 



0<n<N 



Gathering (15. 7p and Lemmas 15.11 and 15.61 we deduce (by modifying if necessary the con- 
stant C from line to line) 



EE E( fc + 1 ) P M su p \Rn\>\h i '\K N = L} 



k>o KeT LeT» °- n ~ N 



<C(l + At) N J2J2 X> + 1) p Ql{ sup \Zn\>jh 1 / 2 ,T N = K} 

Kn<N 4 



k>0 LeT^ KeT 



<C(l + At) N y^$iT k N }(k + iy supQ L { sup IH^I > ^ 1/2 } 

t^ L l<n<N 4 

<^/ 2 exp(CATAt)^(A; + l)^[exp(-^)+exp(-^A_)]. 

fc>0 

By ( 15. 4p and (J53J), we deduce (up to a new value of C) 

< Ch {p+l)/2 NAt{l + NAt) p 

+ C/^ 2 exp(CiVAt) J> + l) d+ v [exp(- J^) + exp(-^L)] . 

fc>0 
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Keeping the inequality h < 1 in mind and comparing the left-hand side below to an integral, 
there exists a positive constant C p , only depending on p and on the same parameters as C, 
such that 



fc>0 

Finally, 



J> + l)*»[„p(- * ) + exp(- * )] < C(l + NM)^'\ 



Theorem 5.7. Assume (14. ip . (15. ip . /i < 1 and At < 9h. Then, for any p > 1, there exists 
a constant C p > 0, only depending on ||a||oo> a > d, K ! V> @ and p, such that, for all N > 1, 



N-l 



h d/2 E k [\X N -X + AtJ2 a ( x n) IT < C p h p/2 exp{C p NAt). 



KeT n=0 



5.5. Analysis of the Numerical Scheme. We now prove the main results of the paper. 

Proposition 5.8. In addition to the assumptions of Theorem \5. % assume that u° is k- 
Lipschitz continuous. Then, for any p > 1, there exists a constant C p > 0, only depending 
on Halloo, a, d, k, r], 9 and p, such that, for all N > 1, 

sup [h~ d ' 2 W u k ~ u(NAt, -)IIW)] VP < C p h 1 ' 2 exp(C p NAt). 

x&Rd K:\e K -x\<hV* 

Proof. By translation, it is sufficient to prove the bound for x — 0. To simplify the 
notations, we set for any random variable Y with values in lR d : 



\Y\u=[hV 2 J2M\ Y \ p ]} 

K&To 



(Above, To is the set of cells K such that |ejf| < h 1 / 2 .) Of course, || ■ \\ p ,t is a norm on 
the space of M. d - valued random variables with a finite moment of order p under every P^, 
K G T . By Theorem 15.71 and by the properties h < 1 and At < 9h, we have 

JV-l 

VA^> 1, \\X N -Z( y NAt,X ) + AtJ2[ a ( X n)-a(Z(nAt,X ))]\\ pTo < C p h 1/2 exp(C p NAt) , 

n=0 

up to a new value of C p . Following the proof of Proposition 13.41 Gronwall's lemma yields 
(5.8) ViV > 1, \\X N - Z(A^At,X )|| pTo < C p h 1/2 exp(C p NAt), 



again for a new value of the constant C p . Following the proof of Theorem 14. 9[ 

N „.(ATA+ WW ~\ 1 / p 

|| t 
KeT a 

for some C > 0. This completes the proof. □ 



[h- d/2 J2 W u k ~ u (NAt, -Wl^k)] 11 " < 4 X n - Z ( N ^ X *)\\ P)% + Cexp(CNAt)h, 



As a by-product, we obtain the L°° estimate announced in Introduction: 
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Theorem 5.9. In addition to the assumptions of Theorem \5. % assume thatu is K-Lipschitz 
continuous. Then, for any p > 1, there exists a constant C p > 0, only depending on Halloo, 
a, d, k, r], 6 and p, such that, for all N > 1, 

sup sup \v% - u(NAt,x)\ < C p /i (1 - 1/p)/2 exp(C p iVAt). 

KeTxeK 



Proof. Consider a given cell K. By Proposition E 



h d/2p inf \u N K - u(NAt,y)\ < \(3h- d l 2 \\u% - u(NAt, OilWj ^ C p h 1 ' 2 ex V {C p N At). 

yeK ^ ' 

Since u(NAt, •) is Lipschitz continuous with Kexp(K,NAt) as Lipschitz constant, we complete 
the proof. □ 
By a regularization argument, we manage to weaken the required assumption on u° in 
Proposition 15. 8t 

Theorem 5.10. In addition to the assumptions of Theorem 5.1, assume that u° belongs to 
BV(M. d ). Then, there exists a constant C > 0, only depending on Halloo, a, d, k, r\, 9 and 
the BV semi-norm ofu°, such that, for all N > 1, 

\\uK-u(NAt,-)\\ L x (K) <Ch 1/2 exp(CNAt). 

KeT 

Proof. The parameters h and N are fixed for the whole proof (with h < 1). The constants 
"C" and "C" appearing below may depend on the BV semi-norm of u°. By [23| Chapter 
5], we know that the semi-norm BV, denoted by || • ||By(Rd\, decreases by convolution. In 
particular, we can find a smooth function u\ : M. d — > R, such that \\u° — tt/lll.L 1 (K d ) — ftl//2 an d 



|| Vn^Hii^d) < ||w°||w(]R d )- We then set, for all x G M. d , u°(x) = \B(0, h 1/2 )\ 1 J B{0A1/2) u° h {x- 
y)dy. We let the reader check that 

(5.9) - M || L i( R d) < h l/2 {\ + \\u°\\ BV(Rd) ), HVm !^!^) < ||it°|Uv(Rd). 

The reason why we introduce u° is the following: the gradient of u° may be locally bounded, 
in the L°° sense, by local L l norms of Vu°. Indeed, for any i6l d and R > 0, 

(5.10) sup |V«°(y)| < |S(0,/i 1 / a )|- 1 ||VttJ|| £ i (flfeBfh i/ a)) . 

yeB(x,R) 

We denote by u(t,x) = u°(Z(t,x)) the solution, at (t,x), of the transport problem with 
u° as initial condition and by = K K (u°(K N )) the corresponding approximate solution at 
time step iV in cell K. It is clear that there exists a constant C > such that 

(5.11) Vt > 0, \\u(t, ■) - u(t, -)\\ L i m < exp(Ct)||^° - u°|Ui 
Similarly, with Proposition 15.31 at hand, 

E \k\\«k-*k\ < E i^i E^K 



KeT KeT LeT 

:_A+\\^~ l„,0 „-0|Wl i nA+\X^\T\\».0 „- 



= E i^ 1 + ^ A *) E tmK - fi £i ^ C 1 + CAt ) E i L n< - u °l 

LeT KeT LeT 

Iterating the procedure, we obtain 

(5.12) Yl \ R W u k - u k\ < (1 + CAtf \ K W u °k - u° K \ < (1 + CAt) N \\u° - u°\\ L i 
KeT KeT 
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By (E2D, (Oil) and (15321) . it is sufficient to investigate E^ e r W^k ~ u(NAt, -)\\li(k)- For 
a cell i^, we obtain by (15.101) and by stability of the solutions to (11.21) . for all y G K, 

\u(NAt,y) -u(NAt,e K )\ = \u°(Z(NAt,y))-u°(Z(NAt,e K ))\ 

< Cexp(CNAt)h sup |Vu°(z)| 

\z-Z(NAt,e K )\<C exp(CNAt)h 



<CeMCNAt)h l - d ' 2 \\Vul\\ LKB{z{NM 



e K ),C*exp(CAfAi)/i 1 /2)), 



for some C > (which may vary from line to line). Integrating with respect to y, we have 
by inversion of the flow Z(NAt, ■) 



£ ||«(iVAt,.)-«(iVAt, e ^)|| Ll(K) 
r 

Cexp(CWAt)/i 1+d / 2 / [|VitX(z)| V l {M (Mt, ejf )|<c«p(CTAtw]^ 



< 



< C exp{C'NAt)h\\Vu° h \\ L i 



It is thus sufficient to analyse Y2k&t \K\\ u k ~ u(NAt, e^)|. In what follows, we fix a point 
i £ K d and we consider a cell K e T x , the set of cells L such that |ex — x| < /i 1//2 . The 
triangular inequality yields (remind that X = ex under Fx) 

\u K -u(NAt,e K )\ 

(5.13) < J2 E 4\ u °k n - w°(^(A r At,X ))|l {Ml / 2 <|^-^(iVA t) Xo)|<(fc + i)ftV 2} ]. 

fc>0 

We are to bound, under F K , the difference — w°(Z(iVAi, X ))| on the set {kh x l 2 < \X N — 
Z(NAt, X )\ < (k + l)h 1/2 } by the gradient of u° and by (k + 2)h 1/2 . (On this set, every point 
in K N is at distance less than (k + l)h 1/2 + h < (k + 2)h 1/2 from Z(NAt,X ).) Since X = e K 
under F K , we know by stability of the solutions to ([U2D that \Z(NAt,X ) - Z(NAt,x)\ < 
Cexp(CA^At)/i 1 / 2 for some constant C > 0. Therefore, up to a modification of C, the ball 
B(Z(NAt,X ), (k + S)h l l 2 ) is included in the ball B k x = B(Z(NAt,x),Cexp(CNAt)(k + 
l)h 1 / 2 ). By (jOOj) . we obtain (under P K ) 

(5.14) sup |Vu°(z)| < C/i-^llV^lUx^). 

|^-Z(7VAt,X )|<(fc+2)/i 1 /2 

By (157T31 and jOg) , 

£ |#||u£-u(jVA*, ejr )| 

< C^ 1+d )/ 2 ^[(A: + 2)||V<|| L1(B , ) £ P K {|X w -Z(iVAt,Xo)| > fc/i 1 / 2 } 

fc>0 L 

Recall the Markov inequality: for a nonnegative random variable Y and two reals a,p > 0, 
P^{F > a} < a' p E K [Y p ]. Choosing Y = h~ l / 2 \X N - Z(NAt,X )\, a = k and p = d + 3, 

31 



and referring to the proof of Proposition 15.81 (see (15.81) ). we obtain for k > 1: 

J2 W K {\X N -Z(NAt,X )\ > kh 1/2 } < k- {d+3) ®K[(hr 1/2 \X N - Z(NAt,X )\) d+3 ] 

KeT x KeT x 

= k~^h- d / 2 \\h- 1 ' 2 (X N - Z(NAt,X )) \\ d d Xl Tx 
< C(k + l)- {d+3) h- d/2 exp(CNAt). 
The modification of k into k + 1 in the last line permits to recover the case k — 0. Finally, 
J2 \K\\u% ~ u(NAt, e K )\ < Ch 1 / 2 exp(CWAt) ]T [(k + 2)(k + \\Vu° h \\ LHBi) ] . 



KeT x k>0 



When integrating the left-hand side with respect to x over M , we obtain by Fubini's Theorem 
a term equal to a constant times h d / 2 ^2 K( - T \K\\u^ — u(NAt, e#)|. When integrating the 
right-hand side, we have by inversion of the flow Z(NAt, 



/ ||Vtt°||ii(Bfc)(ia; = / / \Vu° h \{z)dzdx 

JR d JVL d i_B(Z(AfAt,x),Cexp(CAfAt)(fc+l)?i 1 /2) 

<C'exj>(C'NAt) I [ \Vu Q h \{z)dzdx 

JR d J B(x,C expfCNA^fk+Uh 1 / 2 ) 



' B(x,C exp(C*ArAt)(fc+l)/i 1 /2) 

< C'(k + l) d h d ^ 2 ex P (C'NAt)\\Vu° h \\ L1{Rd) . 
This completes the proof. □ 

6. Annex 

We now show the concentration inequality used to prove Lemmas 15.21 and 15. 61 

Proposition 6.1. Let (O, A, P) be a probability space and d G N \ {0}. Then, there exists a 
constant c > 0, only depending on d, such that, for any M. d -valued martingale (Y n ) n > with 
respect to a given filtration (TC n ) n >o, satisfying Y — and, for all n > 0, \Y n+ i — Y n \ < 1 
and E[|y^ + i — Y n \ 2 \H in \ < v for a deterministic real v > 0, the following holds for all n > 0: 

2 

Vm > 0, P{ sup \Y k \ >u\ <c[exp( ) + exp( )1. 

o<fe<„ J L v cnv' c /J 

Proof. Without loss of generality, we can assume that d — 1. Indeed, 

d 

Vw > 0, P{ sup \Y k \ > u) < VP{ sup \(Y k )i\ > ud- 1/2 }. 

0<k<n ' ~^ 0<k<n 

In the one-dimensional case, it is sufficient to investigate P{sup 0<fc<n Y k > u}, the lower 
bound following from an obvious change of sign. We then apply [7J Proposition 1.6] (choose 
X k = for k > n, a = u and b = nv in the statement of [7, Proposition 1.6]) 

u 2 

P{ sup Y k >u}< exp(-— -). 

o<fc<n 2{u + nv)' 

There are two cases: either u < nv or u > nv. We obtain 

2 

P{ sup Y k >u}< exp(-— -) +exp(--). 

0<fc<n 4nf 4 
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This completes the proof. □ 
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